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Abstract 



This work addresses the question of precisely what features of few body 
models of halo nuclei are probed by elastic scattering on protons at high 
Q\ ', centre-of-mass energies. Our treatment is based on a multiple scattering ex- 

pansion of the proton-projectile transition amplitude in a form which is well 
adapted to the weakly bound cluster picture of halo nuclei. In the specific 
case of n Li scattering from protons at 800 MeV/u we show that because core 
recoil effects are significant, scattering crosssections can not, in general, be 
deduced from knowledge of the total matter density alone. 

We advocate that the optical potential concept for the scattering of halo 
nuclei on protons should be avoided and that the multiple scattering series 
for the full transition amplitude should be used instead. 
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I. INTRODUCTION 



Models of light halo nuclei have been developed ||T]-§] where the few body degrees of 
freedom of a system of loosely bound valence nucleons orbiting around a relatively tightly 
bound core are properly taken into account. Short range, center of mass and some Pauli 
principle effects are often included in these models. 

In this work we develop a multiple scattering expansion of the nucleon-projectile tran- 
sition amplitude for proton scattering from a few body system. When the projectile is 
composed of weakly bound sub-systems a multiple scattering expansion of the transition 
amplitude in terms of 2-body t-matrices describing proton scattering from the projectile 
sub-systems is expected to converge rapidly [|[] . The elastic scattering observables may then 
be derived directly from this expansion. We contrast this with our earlier work []5]|| which is 
based on a multiple scattering expansion of the optical model operator and therefore treats 
the ground and excited states of the projectile on a different footing. The present approach 
is more appropriate for few-body projectiles at high projectile energy. 

Our aim in this work is to understand the nuclear structure features that should be 
incorporated into the reaction mechanism in order to describe elastic scattering of halo 
nuclei from stable nuclei. In particular it is of considerable interest to examine how far 
elastic scattering observables probe correlation effects among projectile nucleons 0||. 

II. MULTIPLE SCATTERING EXPANSION 

We consider the transition amplitude, T, for scattering of a proton from a many body- 
system composed of a small number of sub-systems. We have in mind, for example, n Li 
assumed to be well described by two valence loosely bound nucleons orbiting around a 9 Li 
core. T can be written as a multiple scattering expansion in the transition amplitudes tj 
for proton scattering from each projectile sub-system I |4]]. We ignore explicit reference to 
excitations of the sub-systems, although each t% may implicitly contain effects due to such 
excitations and will certainly do so if, as we shall assume, they describe elastic proton-sub- 
system scattering. In other words our model assumes that we only need refer explicitly to 
excitations of the projectile which involve changes in the relative motion of the sub-systems 
in the projectile. This is consistent with standard few-body treatments of reactions involving 
halo nuclei [p],p~0| . 



The multiple scattering expansion can be written 

T = E^ + E^G £^ + --- (i) 

X X JjX 

where the proton - X subsystem transition amplitude satisfies 

ii = vi + v I G ix . (2) 



The propagator Go contains the kinetic energy operators of the projectile and all the target 
subsystems, Go = {E + — K) . Here E is the kinetic energy, E = ^ — '- in the overall center 
of mass frame, and fi^A is the proton-projectile reduced mass. We ignore the interaction 
between projectile sub-systems in Go (impulse approximation). We note that in the multiple 



scattering expansion eq.(|l]) both elastic and inelastic excitations of the relative motion of the 
subsystems in intermediate states are taken into account. For proton scattering from halo 
nuclei the inelastic channels associated with breakup of the halo nucleus into its sub-systems 
are expected to contribute significantly to the transition amplitude. 

In this paper we truncate the series in eq.([L]) at the double scattering terms. We have 
not evaluated third order terms and we do not claim that they are negligible. They could be 
handled using the techniques of, for example, reference |TTJ. Our purpose here is to assess 
the applicability of the standard approach for proton scattering on light nuclei. We will 
show that inadequacies show up even at the second order level. 

A second aim of our work is to understand the role of various types of correlations in 
elastic scattering from halo systems. In this paper we make a numerical study of the case of 
proton scattering from a n Li projectile at intermediate energies. Our formalism could also 
be applied to p- 6 He scattering which has been studied extensively elsewhere using methods 
which do not use a truncated multiple scattering expansion |T2|] but do not lend themselves 



well to delineating the role of correlations in an explicit way. 

We assume that the projectile wave function can be written as the product of the core 
internal wave function (p c and the wave function of the two body valence system relative to 
the core (p nn (r, R), where r = f 2 — r 3 is the relative position of the two valence bodies 2 and 
3, and R is the vector from the core centre of mass (particle 4) to the centre of mass of the 
valence pair. 

For projectile energies in the intermediate energy region the relative momentum between 
each subsystem pair is small in comparison with the projectile momentum and will be 
neglected wherever it appears. The elastic transition amplitude to second order in the 
proton-subsystem transition amplitudes, involves single scattering terms where the projectile 
scatters from each target subsystem and double scattering terms where the proton scatters 
from one subsystem and rescatters from another. 

A. Single scattering 

The contribution to the single scattering term from proton scattering from one of the 
valence particles, for example particle 2, is given by 

(fc/$|< 12 |fc<$) = (kf<Pnn\il2\ki<Prm) = ^12^12, A)p„(A) (3) 

where p v (A) is defined in terms of the 2-body halo density 

P2 (A 1 ,A 2 ) = jdQ 1 dQ 2V * nn (Q 1 ,Q 2 ) Vnn (Q 1 + A 1 ,Q 2 + A 2 ) , (4) 



by 



p v (A) = p 2 (^-A,^A) , (5) 



where M 23 = m 2 + m 3 ,M 2U = m 2 + m 3 + m 4 , etc. In eq.fl|) (pnniQhQz) is the Fourier 
transform of wave function of the two body valence system relative to the core (p nn (r, R). In 
the case m 2 = 777.3 — "J„ the quantity p v (A) is just the Fourier transform of the probability 



density p(x) of finding a valence neutron at a distance x from the center of mass of the 
projectile as defined by Zhukov et al 0. 

The energy parameter uj\ 2 in eq.(|]) is given by 

miM 34 
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(6) 



and reduces to uj\ 2 = E/2 in the limit of m 4 ^> m 3 , m 2 . 

The contribution to the single scattering term from proton scattering from the core is 



(fc/$|£l 4 |fci$) = (^corc|tl4(^14,A)|^ core )p 2 (0, ^BLA) 



(7) 



where p 2 is defined in eq.(^) and the arguments in eq.(|7p mean that what is involved is the 
density distribution for the motion of the core center of mass, as defined by Zhukov et al [|TJ , 

p 2 (0,#^A) = [ dQ 1 dQ 2 tf* nn (QuQ2)<Pnn(Ql,Q2 + -^ &) • (») 
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and the energy parameter cuu is given by 

(9) 

In the limit of m 4 > 1, u u = E and p 2 (0, ^^A) --> p 2 (0, 0) = 1 so that eq.(^) reduces to 
the expected expression for the proton scattering from subsystem 4. 

Within our model, there are two contributions to the single scattering term. Firstly a 
valence contribution given by the product of the projectile valence system transition ampli- 
tude and p v (A). Secondly a core contribution in which the nucleon-core transition amplitude 
is modulated by the form factor p 2 (0, -^-A) whose departure from unity arises from the 
motion of the core centre of mass about the projectile center of mass. This modulation dif- 
fers from standard applications of the multiple scattering expansion of the optical potential 
operator |6||| that modulate the core matter density distribution pc by that form factor. 

The relevant halo structure information for the single scattering term is thus contained 
in the matter density form factors p v (A) and p 2 (0, -^23_A). 



B. Double scattering 

We next evaluate the double scattering term in the n Li case. We distinguish the terms 
where the proton scatters from the valence neutrons 2 and 3 and the term where the proton 
scatters once from the core and once from a valence particle. In the former case we find that 
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where p 2 (Ai, A 2 ) is defined in eq.([|) and 



Go(q) 



,Ml(23) 



(m 3 kf + m 2 ki 



M- 



+ q\ +ie 



23 



:n: 



In the case of a heavy core, 

}^ (%(?> -^tA) = J dQldQ? Vnn (Ql, Q2) ^Pnn (£1 + ?, Q 2 + A) 

df 2 df 3 e !( S^>\ l( t A - 9 >^ m (f 2) f 3 )| 2 , 
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where ^ nn (f^,f^) = <p nn (f,B) is the valence wavefunction expressed in terms of r 2 and r 3 , 
the position vectors of the 2 valence particles relative to the core. Therefore, this density 
function involves two-body correlations among the valence particles even in the heavy core 
limit. 

The valence system - core double scattering term is given by 



($|ti 2 Git 14 |$) = Jdqt 12 (oJi 2 , ^A + q) {tp com \hl (wu, i^A - q) |^ corc ) 
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where 
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In contrast to the double scattering term arising from the two valence particles, the particular 
elements of p 2 which enters in the heavy core limit is just the one body density of subsystem 
2 in the halo. 



C. Numerical results for Li scattering at 800 MeV/u 

In order to obtain some quantitative idea of the various terms we have identified, we 
have evaluated the multiple scattering expansion for the specific case of proton scattering at 
800Mev/u. from a 3-body model of n Li. For the purposes of the estimate, only the central 
components of the transition amplitudes were taken into account. Coulomb interaction 
effects were not included. 

For the description of n Li we take the Fadeev wave functions of Thompson and Zhukov 
referred in that work as the P3 model. In describing the 9 Li ground state matter density 
distribution we consider a simplified structure model of a Gaussian distribution with a range 
chosen to reproduce the rms radius ||. The first and second order terms were evaluated 
using a NN transition amplitude derived from the Paris potential [0,0 evaluated at the 
appropriate fixed energy parameter with finite mass effects properly taken into account. The 
transition amplitude for proton scattering from 9 Li was generated by an optical potential 
calculated in the single scattering approximation appropriate for intermediate energy elastic 
scattering ||. 



In the evaluation of the second order terms, the propagators were evaluated using the 
eikonal approximation and the principal value term was neglected. For example we use 

GM = -^{ T ^—\ ■ (16) 

An explicit evaluation using gaussian functional forms for the transition matrices and densi- 
ties involved shows that for small scattering angles the ratio of the principal value and delta 
function terms in eq.(|T6|) is less than 1/kiR, where R is a measure of the halo size. This 
ratio is very small in the cases we consider. 

In Fig.l we show the differential cross section for n Li scattering from a proton target 
at 800MeV/nucleon for centre-of-mass in the range we expect to be covered by experi- 
ments (e.g. |§). The dashed curve was evaluated from the single scattering contributions 
eqs.(^),(0). The solid curve includes in addition the double scattering contributions valence- 
valence eg. flTOl) and valence-core eq. ([13]) . We emphasise that in the present context "double 
scattering" means 2nd order in the proton-subsystem t matrix. Terms of all orders in the p- 
subsystem potentials are included in our 1st order terms. The other curves in the figure are 
obtained by taking p2{0, t^-A) = 1 f° r a ^ ^ * n ec L-(0)- This limit corresponds to ignoring 
the relative motion of the core and projectile centres of mass. The dotted-dashed and dotted 
curves correspond to single and double scattering calculated cross section respectively, and 
clearly shows that the inclusion of the relative motion of the core and projectile centres of 
mass has a significant effect in the calculated differential cross section. 

III. DISCUSSION 

In the context of nucleon scattering from conventional stable heavy nuclei one usually 
associates 2-body correlation effects with the double scattering terms which in the present 
case would mean through the 2-body density p 2 (^i 'm^~ ^ in eq.([Tl]). The contribution from 
this to the second order term are very small here, and only valence-core double scattering 
contributions remain relevant. However, that does not mean that the scattering is sensitive 
only to the projectile density p v {A) of eq.(p|), which we might reasonably call "the halo 
density". The scattering involves the halo wavefunction in several other distinct ways: 
Firstly, through the 2-body density P2(0, -^-A) of eq.(||). In the limit of a infinite massive 

core, lim Tn4 ^ oo p 2 (0, -^-A) = 1, but this is a poor approximation in the cases considered. 

Secondly, the halo wave function is involved through the 2-body density P2(jj 2 -q + i^p-A, q) 

in eq. flT3|) . In the limit of a infinite massive core, lim m4 _ >00 p 2 (^ 2 -'?+ jj^-A, q) = p v (A), and 
this limit was found to be a good approximation here. 

There are several consequences which flow from our analysis. 

We have shown that core recoil effects are important. The same claim has been made 
by others but within the framework of formalisms which differ from ours. [|l(| corrects the 
projectile matter density as a whole for recoil effects. One of our points is that we find no 
justification for describing the scattering of protons from a light system such as n Li in terms 
of an optical potential expressed as a nucleon-nucleon transition amplitude, tNN, and a total 
matter density given by the sum of the valence density and a core density modulated by a 
centre-of-mass factor ^2(0, -j^-A). In the first place even in the 1st order term it is the N-core 



transition amplitude which is modulated in this way. Secondly centre-of-mass corrections 
to the second order terms do not have the structure that would arise from iterating the 1st 
order term as would be expected in a "tp" type optical model theory. Eqs. (0),(0) an d ([13]) 
can be made to have this structure if the following 3 assumptions are made: 

(i) ti4 is approximated by its "tp" limit. 

(ii) The average t^N matrices for the core and halo nucleons are assumed equal. 

(hi) The limit m^ — > oo for the core mass is assumed. 

In our calculations we can find no justification for (i) and the inadequacy of (ii) was 
shown very clearly in |5|.|TT|]. We have shown here that (iii) is a poor approximation in 
eq.([8p. In |T7| core recoil effects are taken into account in a way which is consistent with a 
few-body model of the reaction and without making a multiple scattering expansion. It is, 
however, not as transparent as in our formalism how the halo density functions contribute 
to the scattering. An advantage of our approach is that reaction mechanism and structure 
effects are clearly delineated. 

IV. CONCLUSIONS 

We have seen in this work that 2-body correlation effects associated with the double scat- 
tering term are small in the case of n Li scattering from a proton target at 800MeV/nucleon. 
The density distribution of the core centre of mass, P2(0, -^-A), does, however, have a large 
effect on the calculated crosssection. 

We have shown that the halo structure information associated with p v (A) and 
P2(0, -^21. A) does not contribute to the scattering simply combined as a total matter density. 
Thus, a proper treatment of the reaction mechanism for halo nuclei elastic scattering needs 
necessarily to incorporate structure features that go beyond knowledge of the total halo 
matter density distribution alone. 

In summary, we advocate that in microscopic theories of proton scattering from light 
nuclei such as halo nuclei, at intermediate and high energies the multiple scattering series 
for the full transition amplitude should be used and that the optical potential concept should 
be avoided. 
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FIG. 1. Differential cross section for proton scattering from n Li at 800MeV/nucleon. The 
dashed curve was evaluated from the single scattering contributions. The solid curve includes 
in addition the double scattering contributions. The other curves in the figure are obtained by 
ignoring the relative motion of the core and projectile centres of mass in the single scattering term 
and the dotted-dashed and dotted curves are cross sections calculated without and with double 
scattering terms respectively. 



